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^ , ferromagnetic Kondo lattice model on a triangular 
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^ I Abstract. Effect of the coupling to mobile carriers on the 120° antiferromagnetic state is 

^ . investigated in a ferromagnetic Kondo lattice model on a frustrated triangular lattice. Using 

C/3 ' a variational calculation for various spin orderings up to a four-site unit cell, we identify the 

ground-state phase diagram with focusing on the lightly-doped region. We find that an electron 
doping from the band bottom immediately destabilizes a 120° coplanar antiferromagnetic 
order and induces a noncoplanar three-sublattice ordering accompanied by an intervening 
phase separation. This noncoplanar phase has an umbrella-type spin configuration with a 
net magnetic moment and a finite spin scalar chirality. This spin-canting state emerges in 
competition between the antiferromagnetic superexchange interaction and the ferromagnetic 
double-exchange interaction under geometrical frustration. In contrast, a hole doping from the 
band top retains the 120°-ordered state up to a finite doping concentration and does not lead 
to a noncolpanar ordering. 
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, 1. Introduction 

Coupling between charge and spin degrees of freedom is a fertile ground for various fascinating 
. phenomena in strongly correlated electron systems. One of the fundamental models for 

describing the interplay is the ferromagnetic Kondo lattice model, or equivalently, the double- 
exchange (DE) model [1]. The model includes the ferromagnetic Hund's-rule coupling between 
itinerant electrons and localized spins as well as the antiferromagnetic (AF) superexchange (SE) 
^ I interaction between localized spins. In addition to the early researches [HO [3], the rediscovery 

5h ■ of the colossal magnetoresistance (CMR) effect has stimulated intensive studies and many 

interesting properties have been revealed [H [5l [6]. Particularly, in the lightly-doped regions, 
intricate behaviors show up because of the keen competition between the AF SE interaction and 
the ferromagnetic DE interaction induced by mobile carriers. In the early stage of the study, de 
Gennes predicted a spin canting state with a spin-flop type ordering, which smoothly connects 
collinear AF ordering in the undoped insulator and ferromagnetism in a doped metal [3]. Later, 
the scenario has been revisited; it was argued that a phase separation (PS) takes place to hinder 
the canting state and play an important role in the CMR effect [6l [7j . 

Recently, increasing interest has been devoted to the effect of geometrical frustration of the 
lattice structure on the charge-spin coupled systems. This has been inspired by the discovery 
of peculiar transport properties, such as an unconventional anomalous Hall effect, in some 
pyrochlore oxides [8l[9l[T0]. Theoretically, it was shown that, in the ferromagnetic Kondo lattice 



model on a two-dimensional kagome lattice, the coupling to a noncoplanar spin ordering leads 
to the anomalous Hall effect via the so-called Berry phase from a finite spin scalar chirality . 
The idea was extended to other frustrated lattices, such as a face-centered-cubic lattice |12j . 
a triangular lattice [13\ I14j. and a pyrochlore lattice \15\ I16j. Among them, in the triangular- 
lattice systems that we focus on in the present study, it was pointed out that a perfect nesting 
of the Fermi surface at 3/4 electron filling might lead to a noncoplanar four-sublattice ordering 
and the anomalous Hall effect jl3j . The authors performed careful energy comparison among 
possible spin orderings, and found that a noncoplanar four-sublattice spin ordering with a finite 
spin scalar chirality emerges near 1/4 filling in a wider parameter range than in the 3/4 filling 
case [13]. The 1/4 filling phase was later confirmed by Monte Carlo simulation |17l I18| . Such 
triangular-lattice charge-spin coupled systems have also drawn attention experimentally, e.g., in 
some delafossite oxides [1^1 [2U1 [2T| [22] . 

In this contribution, we extend the analysis of the ferromagnetic Kondo lattice model on 
the triangular lattice with focusing on the lightly-doped region. As mentioned above, the 
geometrical frustration brings about new aspects in the charge-spin coupled systems, but the 
keen competition between the AF SE and the ferromagnetic DE interactions in the lightly-doped 
region has not been studied in detail. On the triangular lattice, the AF SE interaction stabilizes 
a three-sublattice 120° AF order when the electron band is empty or fulfilled. We clarify the 
effect of both electron and hole doping to the 120° AF states by a variational calculation of the 
ground-state energy for various spin states up to four-sublattice ordering. In the electron doped 
case, we find that the 120° order is immediately collapsed by doping and a three-sublattice 
spin canting phase appears with accompanying PS. In contrast to the unfrustrated cases, this 
canting state has a noncoplanar umbrella-type spin configuration, which shows a finite scalar 
spin chirality together with a net magnetic moment. On the other hand, in the hole doped case, 
the 120° AF state persists up to a finite doping rate and PS takes place to a ferromagnetic metal; 
a noncoplanar canting phase is not observed. The contrastive behavior between the electron and 
hole doping is discussed from the change of the density of states for the spin ordering patterns. 

2. Model and method 

We consider the ferromagnetic Kondo lattice model on the triangular lattice. The Hamiltonian 
is given by 

'^ = -t J2 + h.C.) -JrJ2 4,a^»l3Ci,(3 ■ + JkJ2 Si ■ Sj, (1) 

(i,j),a i,a,^ (i,j) 

where c|„ {ci,a) is a creation (annihilation) operator for a conduction electron with spin a at 

site I, ffap = {(^ap^'^ais^^afs) ^ vector of Pauli matrices, and Si is a localized spin on site 
i. Here, t is the transfer integral between the nearest-neighbor sites on the triangular lattice, 
Jh is the Hund's-rule coupling between conduction electrons and localized spins, and Jk is the 
antiferromagnetic (AF) superexchange (SE) interaction between localized spins. We take t = 1 
as an energy unit and consider classical spins for Si with \Si\ = 1. 

We investigate the ground state of the model given by Eq. (1) following the variational 
method taken in Ref. [T^. Namely, we compare the ground-state energies for different ordered 
states of the localized spins and determine the most stable ordering, while varying the electron 
density n = ^ J^iai'^l a'^i,a) is the total number of sites), Jh, and Jk- In the calculation, we 
consider 13 different types of ordered states, up to a four-sublattice unit cell, as shown in Fig. 1 
in Ref. |14j . For the states with spin canting, we optimize the canting angle 9. 
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Figure 1. Ground-state phase 
diagram as functions of n and Jh 
at Jk = 0.03. Labels for the 
phases are common to those in 
Ref. [H]. (3b) is the umbrella- type 
noncoplanar phase appearing in 
between the ferromagnetic metallic 
state (la) and the phase separation 
(PS). The thick black line at n = 
indicates that (3a) 120° coplanar 
ordered phase is stabilized at zero 
doping. 



3. Result and Discussion 

Figure [1] shows the result of the ground-state phase diagram at Jk = 0.03 in the lightly-doped 
region near n = 0. For comparison, the results at smaller Jk are found in Fig. 3 in Ref. ^14j . 
A new feature emerges in the low-density region in Fig. [H which was not seen in the previous 
results for smaller Jk. That is, a spin-canting phase (3b) appears on the lower-n side of the 
ferromagnetic metallic phase (la). This state has a three-sublattice ordering with an umbrella- 
type noncoplanar spin configuration, as schematically shown in the figure [see also the inset 
of Fig. [2]^a)]. Spin canting itself is expected in general in lightly-doped region because of the 
competition between the AF SE interaction Jk and the ferromagnetic DE interaction [3]. For 
example, in the cubic lattice plausible canting state, which naturally connects two- 

sublattice collinear AF state and doped ferromagnetic state, is a simple coplanar spin-flop state. 
However, it is not trivial what type of spin-canting order is realized in geometrically-frustrated 
systems. In the present triangular-lattice case, there are several possibilities of different spin 
canting orders, as listed in Fig. 1 in Ref. [H]; even for the three-sublattice orders, there are 
three different types of canting states, (3b) umbrella, (3c) coplanar cant, and (3d) 2:1 cant (see 
Fig. 1 in Ref. [H]). Among the possibilities, the umbrella- type one is energetically favored in 
the present model. Note that the umbrella state has a finite scalar spin chirality in each triangle 
unit, but this spatial configration does not lead to the anomalous Hall effect. 

The umbrella phase is developed by increasing Jk in the region near the border of the 
ferromagnetic phase to PS as Jk increases. PS occurs to the n = state, in which itinerant 
electrons are absent and the localized spins are 120° ordered by Jk- It is noteworthy that 
the umbrella phase is always accompanied by PS. The situation is demonstrated in Fig. [2j 
Figure EJa) shows the energy per site at a fixed canting angle 9 for various values of 9. Here, 
^ = 0° and 9 = 90° correspond to the ferromagnetic phase (la) and the three-sublattice 120°- 
ordered phase (3a), respectively. As 9 decreases from 90°, the energy at n = increases because 
of the loss of the SE energy. At the same time, the initial slope dE/dn\n=o becomes steeper 
because the kinetic energy of electrons becomes larger for smaller 9 via the DE mechanism. 
Hence, the envelope of the energy curves E(n, 9) becomes convex upward in the low-density 
limit. The convex behavior is a sign of PS. The phase-separated region is determined by drawing 
a tangent line from the energy at n = 0, E{n = 0,9 = 90°), to the envelope of the energy curves 
E(n > 0,9 < 90°). In the present case at Jk = 0.03 and Jh = 10, the tangent line touches the 
envelope at n ~ 0.06 with the energy curve for 9 ~ 50°; i.e., the system exhibits PS between the 
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Figure 2. (a) Energy per site of tiie umbrella-type noncoplanar ordered state. Different curves 
represent the energies for a fixed different canting angle 9. The data are calculated at Jk = 0.03 
and Jh = 10. Note that 6 = 0° and 6 = 90° are equivalent to (la) ferromagnetic phase and (3a) 
three-sublattice 120°-ordered phase, respectively. The straight black line shows a tangent line to 
define the phase-separated region. The insets show schematic pictures of the spin configuration 
of the umbrella state and the triangular lattice, (b) n dependence of 9 which gives the lowest 
energy. The gray region denotes the phase-separated region. See the text for details. 



120°-ordered state at {n, 9) = (0, 90°) and the umbrella-type canting state at (n, 9) ~ (0.06, 50°). 
Similar behavior is observed for all the regions of the umbrella phase, which accounts for the 
reason why the umbrella phase always appears with accompanied by PS. 

The canting angle 9 which gives the lowest-energy envelope changes continuously in the 
umbrella phase, as seen in Fig. [5]|^a). Figure [5]|^b) gives a plot of the optimized as a function 
of n. The result shows that, as n increases, 9 starts from about 50° at the border of PS at 
n ~ 0.06, and gradually decreases to zero when approaching the boundary to the ferromagnetic 
phase at n ~ 0.13. 

The value of the optimized 9 depends not only on n but also on Jh (and Jk). The behavior 
is summarized in Fig. Oat Jk = 0.03. As shown in the figure, 9 decreases as n and Jh increase. 
This is reasonable since larger n and Jh enhance the DE ferromagnetic interaction which favors 
smaller canting angle 9 (larger net magnetic moment). This change of 9 clearly demonstrates that 
the noncoplanar umbrella phase is stabilized in the competition between the AF SE interaction 
Jk , which is dominant at less carrier doping and small Jh , and the DE ferromagnetic interaction 
enhanced by doping and Jh. 

In contrast to the above results for the electron doped case, a hole doping from the band 
top n = 1 does not lead to a noncolpanar canting state. In the hole-doped case, 120° coplanar 
ordering persists up to a finite doping concentration, and PS takes place between the doped 120° 
AF metal and a ferromagnetic metal when Jh is sufficiently large. The situation is unchanged 
against the value of Jk; a typical phase diagram for n ~ 1 is found in the result at Jk = 0.01 in 
Fig. 3 in Ref. [H]. 

The reason why the 120° AF order is stable against the hole doping is understood from the 
density of states (DOS). Figure H] shows DOS for the 120° coplanar AF phase, in comparison 
with typical DOS for the umbrella canting states including the ferromagnetic state. The data 
are calculated at Jh = 10 and Jk = 0. The result shows that the band top energy e = 13t 
is unchanged for 9. This is confirmed by analytical calculations. In particular, a consideration 
in the limit of Jh ^ oo provides simple understanding of the form of DOS in terms of the 
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Figure 3. Canting angle 6 in the 
umbrella phase as functions of n and Jh 
at Jk = 0.03. The ferromagnetic phase 
(9 = 0°) and the 120°-ordered phase 
{9 = 90°) are also shown by black and 
white, respectively, whereas other phases 
are omitted (the hatched areas). The 
gray areas are the phase separation. 



renormalization of the transfer integrals and the flux introduced by the Berry phase [21 I23j . 
While the band top energy is unchanged, the value of DOS at the band top is smaller for 9 > 
compared to the 9 = ferromagnetic state. For these situations, the kinetic energy by hole 
doping becomes lower for 9 > states. Among the 9 > states, the 120° AF state has the 
lowest AF SE energy, and therefore, the 120° AF ordering remains stable against the hole doping 
for Jk > 0. 
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Figure 4. The density of states 
(DOS) for the umbrella-type canted 
states including the 120° coplanar AF 
state and the ferromagnetic state. The 
data are calculated at Jh = 10 and 
Jk = 0. 



DOS in Fig. U] also accounts for the reason why electron doping from n = immediately 
collapses the 120° AF state and causes PS to a canted state. As shown in the figure, the 
band-bottom energy is lower for smaller 9, and hence, the electron doping to the 120° AF state 
costs the energy than the canted states. At n = 0, the 120° AF ordering becomes stable by a 
small energy gain from the AF SE coupling, — 3Jk per site. The energy gain is compensated by 
infinitesimally small doping to a canted state with smaller 9 because of the lower band bottom. 
Therefore the electron doping leads to a sudden PS to a canted state, as demonstrated above. 

The contrastive behavior between electron and hole doping is expected to be observed around 
an insulating state at half filling n = 0.5. At n = 0.5, the 120° coplanar AF ordering is stable 
and the system becomes insulating for Jh ^ 2 (see Fig. 3 in Ref. [13 )• In the insulating state 
for larger Jh, DOS just below and above the gap behave similar to the band top and bottom, 
respectively, as shown in Fig. HI In fact, as a preliminary result, we find a similar contrastive 



behavior around n = 0.5; a noncoplanar canting state appears against electron doping with 
PS, while the 120° AF order retains up to a finite hole doping. The details will be reported 
elsewhere. 

Finally, let us discuss briefly the effect of the band structure on these properties by including 
the next-nearest-neighbor hopping t' in the model ([1]). We find that, when t' has an opposite 
sign to t, 6 dependence of DOS at the band edges is qualitatively unchanged. For t' with the 
same sign as t, the band top becomes higher for smaller 6, while the band bottom remains lower 
for smaller 6, at least, for small t'. Hence, we expect the similar results for the models including 
t' irrespective of its sign. Further analyses including the detailed energy comparison among 
different spin states will be reported elsewhere. 

4. Summary 

We have investigated the lightly-doped regions of the ferromagnetic Kondo lattice model on the 
triangular lattice by variational calculations. We have found that a three-sublattice noncoplanar 
phase is stabilized by electron doping from the band bottom. The noncoplanar spin canting phase 
has an umbrella-type spin configuration with a finite spin scalar chirality in each triangular 
plaquette. This phase always appears accompanied by a phase separation to the 120° coplanar 
antiferromagnetic insulating state at zero doping. On the contrary, hole doping from the band 
top does not lead to such a noncoplanar spin canting state. These peculiar properties are 
the consequence of geometrical frustration in the charge-spin coupled system in which the 
antiferromagnetic superexchange interaction and the ferromagnetic double-exchange interaction 
compete with each other. 
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